Emergence of Artificial Photons in an Optical Lattice 
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We establish the theoretical feasibility of direct analog simulation of the compact U{1) lattice 
gauge theories in optical lattices with dipolar bosons. We discuss the realizability of the topological 
Coulomb phase in extended Bose-Hubbard models in several optical lattice geometries. We predict 
the testable signatures of this emergent phase in noise correlation measurements, thus suggesting 
the possible emergence of artificial light in optical lattices. 
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Introduction. Cold atomic gases in optical lattices 
[l[ have provided unprecedented flexibility in designing 
and studying coherent and correlated condensed mat- 
ter systems. To date, experiments on strongly corre- 
lated, bosonic lattice models have realized the on-site 
Hubbard-type interaction U , comparable to a nearest- 
neighbor hopping t. Consequently, Several strongly cor- 
related bosonic phases with short range interaction have 
been studied However, with the recent obser- 

vation of Bose-condensation of Chromium 6] , which has 
a magnetic dipolar interaction, a sizable nearest-neighbor 
interaction V , albeit anisotropic, is now within experi- 
mental reach. Moreover, since U can be tuned using Fes- 
hbach resonances U /V can, in principle, be made to 
vary over a wide range. Furthermore, exciting develop- 
ments [3| in cooling polar molecules offer the possibility 
of strong and tunable electrical dipole moments. 

Recently, it has been proposed [9|, that a four-site 
'ring-exchange' interaction can also be implemented al- 
lowing the simulation of a J7(l)-lattice gauge theory, 
which has various exotic topological phases, among them 
a 3D U{1) Coulomb phase with an emergent massless 
photon mode. In this Letter, we show that by sim- 
ply implementing a Hubbard model with an additional 
strong nearest-neighbor interaction on some special lat- 
tices in both 2D and 3D one can efficiently simulate 
a C/(l)-lattice gauge theory. The dipole interaction, 
Vdd{R) = — 3cos^(0)]/i?'^, where d characterizes the 
dipole moment and R is the inter-dipole separation, is, in 
general, anisotropic and depends on the angle ((> between 
the vectors defining the parallel dipole orientation and 
the bond between the two dipoles. This interaction needs 
to be made isotropic for our purpose. This is easily done 
in 2D by simply aligning the dipoles perpendicular to the 
plane. In 3D, this is possible only on some lattice and we 
take the pyrochlore lattice where this can be achieved. 
Once the gauge theory is simulated, the Coulomb phase 
can be accessed by the appropriate tuning of interaction 
parameters. The existence of a gauge theory automati- 
cally implies that the conventional insulating phases [lol | 



are not the only phases possible. 

We predict and discuss the observable signatures of 
the Coulomb phase, specifically, how the emergent pho- 
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ton mode, giving rise to artificial electrodynamics 
llSll . can be detected in noise correlation measurements 
[l4l . [lil [lit . We note, in passing, that other fraction- 
alized insulators of the gauge theory, where the elemen- 
tary excitations carry fractional boson 'charge' (boson 
'charge' in the present context of neutral bosons means 
boson number), can also be accessed in other regions of 
the parameter space. These, and experimental implica- 
tions thereof, will be discussed in a future work [l7|. 

Dipolar Bosons in Optical Lattices. We study the 
Coulomb insulating phase of dipolar bosons in optical 
lattices, specifically the Kagome (2D corner-sharing tri- 
angles) and pyrochlore (3D corner-sharing tetrahedra) 
lattices, see Fig. [T] Optical lattices formed from the in- 
tensity extrema of standing wave lasers can be used to 
generate a variety of lattices with non-trivial primitive 
unit cells using the superlattice technique [18], originally 
proposed to generate two dimensional lattices. Using 
this technique one can, in principle, generate a well de- 
fined set of intensity maxima defining the Kagome and 
the regular triangular lattices through angular interfer- 
ence of several beams with the same wavelength. Re- 
markably, the potentials defining these two lattices only 
differ by a phase. By generalizing this technique to 
three dimensions one can create a three-dimensional py- 
rochlore lattice by alternately stacking the triangular and 
Kagome lattices. Consider the following intensity pattern 
generated by counter-propagating, red-detuned lasers: 
F,{r) = [cos(fc,-f+$,) + 2cos(fc,-r/3-f «'j/3)]2. With po- 
larization (or frequency) mismatches, intensities can be 
added to generate the following three-dimensional poten- 
tial with maxima defining a pyrochlore lattice: Vp{r) cx 
^2^=1 -^i^ii^ ■ The relative intensities and wavevectors 
are given by: Bi^4 = I, B4 = 2Bi, ki = fc(l, 0,-5/7), 
k2,3 = ^(1,±n/3/2,-1/7), and fcg - ^4 = fc(0, 0,3/7). 
Here we have defined 7 = 4[1 - pV(3p'^)]^^^- The rel- 
ative phases of the standing waves, $^=2-4 = and 
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FIG. 1: In the left figure dipolar bosons sit on the sites of the 
2D Kagome lattice, the black circles labeled 1-6. The hexago- 
nal dual lattice, the lattice formed by the centers of the corner- 
sharing triangles (gray sites), is labeled a — f. In the right 
figure, for simplicity, we show only two corner-sharing tetra- 
hedra of the pyrochlore lattice, with sites shown as spheres. 



<i>i = $5 = — 37r/2, must also be fixed. To generate 
Vp, we require phase matcliing between beams with vary- 
ing spatial periodicities. The fact that the beams must 
also be tuned, energetically, to internal states of con- 
stituent bosons (which normally fixes the periodicity) 
poses an experimental challenge. However, recent experi- 
ments have in fact demonstrated multiple spatial period- 
icities through polarization-selective, angled interference 
and Doppler-sensitive multi-photon transitions [l^, 
separately. A combination of both techniques could of- 
fer a versatile system capable of realizing the standing 
wave arrangement used to produce Vp. Further details 
discussing the parameters and laser geometry for these 
lattices will be discussed elsewhere |17| . 

The interaction between dipolar bosons confined to op- 
tical lattices arises from two terms: the short range, s- 
wave, interaction, Vs{R) — gS'^^^R), and the long range 
dipolar interaction Vdd{R)- In what follows we place the 
system in a uniform external field, pointing along the 
stacking direction (the [111] direction for the pyrochlore 
lattice) to orient the dipoles perpendicular to the basal 
plane. We also assume that the interaction strength g 
is tunable through Feshbach resonances, allowing tun- 
ability of the ratio Vs/Vdd- To make the dipolar interac- 
tion, Vdd{R), isotropic we need a pyrochlore lattice where 
the tetrahedra are shrunk in the direction perpendicu- 
lar to the basal plane, that is, p' < p, where p and p' 
are the tetrahedron side lengths within and out of the 
basal plane, respectively, see Fig.[T] The fact that the 3D 
nearest- neighbor interaction can be made isotropic along 
the bonds by such slight engineering of the tetrahedra is 
crucial, and so the pyrochlore lattice is the geometry of 
choice. 

We now consider a simple model of dipolar bosons 
confined in deep optical lattices in the Hubbard limit. 
The resulting hopping term (tunable through the opti- 
cal lattice depth) is nearest-neighbor while the interac- 
tion has both on-site and extended terms llOl. We omit 



the smaller next-nearest neighbor interaction terms gen- 
erated by the long range part of the dipolar interaction. 
Working in the insulating limit we may take the hopping 
energy gain (in units of the on-site energy) to be much 
smaller than the average number of particles per site. We 
then approximate the system by a quantum rotor model: 

H ^ -t'Y^h\bi, +h.c. + U^nf + V^ninv. (1) 

(IV) I ill') 

Here {W) denotes a nearest-neighbor pair of sites, = 
exp (i9i) is the bosonic creation operator and n; = d/idOi 
counts the excess number of bosons at a site, /. The dom- 
inant contribution to U arises from the s-wave interaction 
while the dipolar interaction supplies the largest contri- 
bution to V . We assume that both U and V are large, 
and, in what follows, we will require U to be slightly 
larger than V . Here and in the following we omit the 
renormalized chemical potential term. 

Mapping Onto Gauge Theory in 2D Kagome lattice. 
Eq. [T] can be written in a slightly different form more 
amenable to mapping onto a gauge theory. We take out a 
piece, 2C/o, of the on-site interaction and group it together 
with V . Defining N/^ = X]/gA ^'^'^ u = U — 2Uo , we 
find, for V = 2Uo, 

H ^-tJ2bX+h.c.+uJ2^f + UoJ2^i- (2) 

(W) I A 

Here, sums over all triangles in the lattice, and A^a 
counts bosons on a given triangle. In this form, we call 
this Hamiltonian a 'cluster' Hubbard model where the 
fluctuations in the number of bosons on each local cluster 
(triangle) costs energy. This property enables a natural 
map to a gauge theory which enforces a local constraint. 

Eq. [5] describes the same Hamiltonian as the one pro- 
posed before [ll|, [l2| for a microscopic model of bosons 
on a square lattice and its 3D generalization, the corner- 
sharing octahedra. We focus here, however, on the 
Kagome lattice, for which the same Hamiltonian has 
an exact implementation in terms of a nearest-neighbor 
Hubbard model for U ^ V . The mapping of Eq. [5] on a 
(2-1- 1)-D U{\) lattice gauge theory is essentially the same 
as that on a square lattice. The details of the mapping, 
however, are different, which we briefly describe below. 

A mapping onto gauge theory becomes possible when 
C/q is the largest energy scale in the problem. The boson 
number on each triangle is then locally conserved. This 
implies Aa = on each triangle, meaning that the fluc- 
tuation in the boson number around the mean value set 
by the chemical potential vanishes. Working in the limit 
t^u and doing degenerate perturbation theory in t 
and M, we find, to first order in u and third order in i, 
the following effective Hamiltonian in the ground state 
sector A'a = 0, 

H,a = uY,n^i-i^Y.^\b2b\hAK (3) 

; ^0 o 
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Here, J2o sums over all hexagons inscribed in the tri- 
angles, and 1 through 6 indicate the corners of one such 
representative hexagon, see Fig. [TJ This effective Hamil- 
tonian is derived by locally conserving the number of 
bosons on each triangle. If a boson hops from site 6 to 
site 5, then another boson must hop from site 4 to site 
3 to conserve the number of bosons on the triangle with 
one side defined by the bond (4 — 5), so that Na remains 
zero on that triangle. In turn, as is clear from Fig. [l] a 
third boson must hop from site 2 to site 1 to preserve the 
same constraint. Finally, since the tt-term of Eq. [5] does 
not fluctuate TVa, it comes in linear order in Eq. [3l 

To see that this Hamiltonian exactly maps onto a lat- 
tice gauge theory on the dual lattice, we consider the dual 
lattice defined by the centers of the triangles. This is a 
hexagonal lattice whose sites we denote by r, r' etc. The 
original Kagome sites now fall on the links of the dual 
lattice. We take the representative hexagon, a through 
/, as shown in the left panel of Fig. [l] and identify the 
bosons on Kagome site 1 with the link (a5) etc. 

On the dual lattice, by a change of notation, HcS reads. 



(rr') 



'-^YpJ2blAA,Keblfbfa- (4) 
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Next, noticing that the hexagonal lattice is bipartite, we 
define a field a^r' — Orr' if r G A and r' G B and a^r' — 
—Orr' li r G B and r' S A, where A and B are the 
two interpenetrating sublattices, and brr' — exp {iOrr') as 
discussed after Eq. [1] We also define the conjugate field 
variable err' — nrr' if r S A and r' € B and err' = —rirr' 
if r e _B and r' £ A {e and a are conjugate fields since 
n and 9 are). Finally, elevating a and e to vector fields 
o-ra = o,r,r+a, ^ra = Br.r+a, whcrc a indicates the nearest 
neighbor vectors, it is straightforward to see that iJeff can 
be written as the Hamiltonian for the (2 -I- 1)-D compact 
U{1) gauge theory [20| on the dual lattice. 
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cos(V X a). 



(5) 



o 



with the Gauss's law constraint Nr ~ A^a = 
X^r'er ^rr' — ?7rV.e, and rjr ~ 1 if r IE A and rjr = —1 if 
r E B. We emphasize that Eqs. are simply rewrit- 

ten forms of Eq. [3] on the dual lattice. Recall that a 
large Uq implied the constraint on A^a in Eq. [5] leading 
to Eq. [3] Hence, for an optical lattice implementation 
of the gauge theory, one simply needs to implement the 
Hamiltonian in Eq.[T]on appropriate lattices with both U 
and V large, with U slightly larger than V. The result- 
ing low energy theory is automatically a gauge theory on 
the dual lattices, and so allows the exotic Coulomb phase 
and other fractionalized phases 111 Ip . [20l . [2lj in addi- 
tion to the conventional insulators [lO|. In 2D, however, 
it is well known that the Coulomb phase is unstable at 
long length scales and is smoothly connected to the con- 
ventional Mott insulator. 



Mapping Onto Gauge Theory in 3D pyrochlore lattice. 
We start with formally the same Hamiltonian as Eq. [U 
but on a pyrochlore lattice. From Eq. [1] one straightfor- 
wardly gets the 'cluster'-Hubbard model, Eq. [21 with N/\ 
replaced by Nt, the number of bosons on a tetrahedron. 
The mapping to gauge theory on the dual diamond lat- 
tice, which is bipartite, but not Bravais, is identical and 
one ends up with Eq. [5] as the final effective theory. To 
have a description of the dual lattice in terms of a Bravais 
lattice, we consider below only one FCC sublattice of the 
diamond lattice (centers of only one class of tetrahedra 
on the pyrochlore). Note that by counting the centers 
of only one class, and allowing an index a to indicate 
all four corners of a tetrahedron from its center, we can 
count all the pyrochlore sites (diamond links). 

In (3 -I- 1)-D, the Coulomb phase, which is an insula- 
tor, is described by the Gaussian expansion of the co- 
sine term of the gauge theory Hamiltonian, and is sta- 
ble for t ^ y'uU§, since the topological defects - the 
monopole configurations of the gauge field - are sup- 
pressed in the Coulomb phase. This is analogous to 
the gaussian expansion of the Hamiltonian of the XY 
model in the Kosterlitz-Thouless phase where the topo- 
logical defects - the vortices - are suppressed. Notice 
that a Gaussian expansion of the cosine term yields a 
Hamiltonian formally the same as that in classical elec- 
trodynamics, with an associated gapless 'photon' mode. 
However, this is a gapless mode of the gauge field of the 
bosonic field operators, and not of an external electro- 
magnetic field. This mode, which emerges in the low en- 
ergy theory, should therefore be observable in the long- 
wavelength boson density-density correlation functions, 
and hence in noise correlations, given by the correlation 
functions of the e-fields (Recall that e^r' is simply re- 
lated to Urr')- This mode distinguishes itself from the 
conventional phonon mode resulting from spontaneously 
broken symmetry in strongly correlated, bosonic models 
because it appears in the absence of spontaneous symme- 
try breaking and is therefore an emergent phenomenon. 
Furthermore, this mode appears in a range of parame- 
ters {Uq ^ t ^ -y/ uUq) intermediate between that for 
a superfluid {t 3> Uo,u) and the conventional insulators 
(J7o,u> 0- 

Noise Correlations in the Coulomb Phase. We now 
consider the long wavelength behavior of noise correla- 
tions in terms of the e-field correlation functions. To 
compute the correlation function of the e-fields, we work 
in the Coulomb gauge where Or = and V • a = 0. We 
then expand the cosine in the second term of Eq. [5] to 
quadratic order and use e = Working in the con- 
tinuum limit, the a- field correlator at momentum k and 
(imaginary) frequency lu is given by. 



{ai{k,uj)aj{-k, -uj)) 



(6) 



where = (6i^)/(wJ7Q ). This reveals the mass-less 
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FIG. 2: Expected behavior of the normahzed noise correlation 
function with in the Mott insulating (dashed curve) and 
Coulomb (solid curve) phases around the first Brillouin-zone 
center of the FCC lattice. The core of the Mott insulating 
peak is suppressed due to the c\q\ non-analiticity provided 
by the gapless 'photon' at long wavelengths in the Coulomb 
phase. The behavior, at large q, is expected to cross over 
to that of the Mott insulator. The plot is for a finite sized 
system with 100 sites along each direction with a small q slope 
c = 40 (with the appropriate form factor set to unity) and an 
arbitrarily chosen cut-off at q^ = Stt/SOO. 



emergent 'photon' mode. The equal-time correlator 
{era&r'a'), where a, a' are unit vectors for a tetrahedron, 
can be found from here. This gives the boson density- 
density correlation function. We find. 



aittj c\k\{Sij 

k 



ik-{r — r ) 



(7) 



A straightforward evaluation of this integral in the 
continuum gives the special angular structure of the 
Coulonib phase density correlation function at large dis- 
tances [2^. However, as described below, for noise cor- 
relation, one is interested in the function in momentum 
space itself. 

In time of flight imaging of an insulating state of spin- 
less particles, averaging the noise between shot-to-shot 
images of the particle distribution released from optical 
lattices reveals a quantity proportional to the following 
second order correlation functionlljll: 



Jq-lf+a-f' -a') 



7^ Ct ^ v' Ot.' 



(8) 



where we, for simplicity, omit the q = Q delta func- 
tion due to normal ordering. In a Mott insulator, 
which should be compared with the Coulomb insula- 
tor, {riranr'a') IS a constant, and so, for an infinite 
system N(q) ^ S«/('7)'^('f ^ ^n) where g„'s are the 
reciprocal lattice vectors of the FCC lattice, and f{q) 
is the form factor of the tetrahedron basis, f{q) = 



Sao' ^^P I*?' ~ '^')]- This produces 5-function peaks 
at the reciprocal lattice vectors which are the centers of 
the Brillouin zones. For the Coulomb phase, we know 
the behavior of {nraUr'a') only for large |r — r*|. At 
smaller length scales, the calculation is quite involved 
since the full tetrahedron lattice correlators have to be 



used [17|. The long distance behavior, however, is suf- 
ficient to give the leading singularities at the Brillouin 
zone centers. From the integrand of Eq. [7] we see that, 
close to the zone centers, the leading singularity is a cusp, 
N{q)^c\q\,see Fig.H] 

In conclusion, we have shown that compact U{1) lat- 
tice gauge theories can be simulated in suitably designed 
cold atom optical lattices using dipolar bosons. Exciting 
exotic phenomena, the Coulomb phase with emergence 
of artificial photons being an example, could be studied 
experimentally following our suggestion. 
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